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SECTION  I 


INTRODUCTION 

This  volume  documents  the  mathematical  methodology  in  the  conceptual 
development  of  the  numerical  and  analytical  relationships  for  the  computa- 
tional trajectory  simulation  models.  These  models  were  exercised  for 
specific  targets/sensors  complex  scenarios  tempered  by  launch  site  locations, 
sensor  visibility,  and  performance  parameters.  These  models  were  represented 
by  a modification  of  a previously  developed  RADC  Trajectory,  Simulation 
Computer  Program.  It  is  to  be  noted  that,  for  the  purpose  of  clarity, 
certain  mathematical  developments  of  these  analytical  expressions  do  not 
employ  the  common  scalar  notation  of  a vector.  However,  this  scalar  form  is 
equivalent  to  the  vector  symbol  without  the  arrow. 
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SECTION  2 


TRAJECTORY  SIMULATION 

The  computer  program  used  for  the  calculation  of  trajectories  is  a 
special  perturbation  computer  program  generated  at  RADC.*  It  was  run  on 
the  RADC  HIS-635  computer.  The  double  precision  program  provides  sensor 
metric  data  and  their  time  rates  for  specific  classes  of  trajectories.  These 
simulated  trajectories  encompassed  reentry  type  vehicles:  ICBM,  SLBM,  classes 

and  the  interceptor/synchronous  target  scenarios.  These  Hohmann  and  non- 
Hohmann  ballistic  flight  vehicles  were  launched  from  specific  geographic 
locations,  i.e.,  at  booster  burnout,  to  impact  in  specific  locations.  The 
interceptor/target  scenarios  included  classes  of  Hohmann  and  non-Hohmann 
trajectories  that  can  be  simulated:  (i)  direct  ascent  maneuvers,  (ii)  low- 

alt  itude-dwell  scenarios  (elliptical  transfer)  to  high-altitude-dwell 
scenarios  (bielliptical  transfer).  The  time  history  of  each  vehicle's  dynamic 
behavior  is  influenced  by  a rotating  spherical  earth  potential,  standard 
atmosphere,  earth  wind  and  a variable  ballistic  coefficient.  A discussion  of 
the  main  subroutine  groupings  used  are: 

( 1 ) ICF’RBM  With  an  input  of  a set  of  initial  boundary  values  unique  to  each  vehicle,  the 
tree-flight  position  and  velocity  vectors,  relative  to  an  inertial  geocentric  coordinate  system,  are 
calculated  via  the  transformation  equations: 

dr/dt  |G  = V 0 cos  7 I 

d0/dt  |0  =(VG  sin  7 cos|3)/r0  J 

dX/dt  |Q  = (VG  sin  7 sin  0)/ro  cos  0 1(11 

where  the  free  flight  reentry  angle  (WRT  local  vertical),  heading  angle  (WRT  North)  and  the 
magnitude  of  the  velocity  vector  (relative  to  burnout)  are  given  respectively  by  7.  (3,  and  V0  (as 
illustrated  in  Figure  2-1 ).  In  addition,  the  salient  feature  of  this  routine  is  to  provide  a corrected 
set  of  initial  boundary  conditions  corresponding  to  a specified  impact  location. 

(2)  EQTNX  With  respect  to  the  inertial  geocentric  coordinate  system,  the  vehicle's 
acceleration  components  are  calculated.  These  components,  in  spherical  coordinates,  are  defined 
by  the  following  set  of  2nd  order  non-linear  differential  equations  of  motion. 

d2r/d2t  — r02  (ruis!  0)  S22  + p/r2  = Oi  I 

d(r20)/dt+(r2  sin  20 ) S 2 2 / 2 = Q2  | 

d(r2  S2  cos20)/dt  =Qj  J (2) 

where  Qj  represents  the  nongravitational  force  arising  from  atmospheric  drag  i.e., 

01  =-pZdVr 

02  =-pZ<jVr20 

(J)  =-pZtjVr2  (X  >!  cos2 0 (3) 

* I Ills,  (I.A.,  et  al  , “Simulation  Program  for  Three  Degrees  of  Freedom  Trajectories,” 

EMA-TM-66-5,  Sep  66. 
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where;  p = 

n = 

'•d  = 

(unit  mass) 

wc  = 

ijj  = 


atmospheric  density 
u;  + X 

l'jA/2  = ballistic  coefficient 


earth’s  angular  velocity 
earth  wind 


North 


1 


r 


Figure  2-1  Vehicle’s  Velocity  Components 

(3)  KK(i  Using  a 4tl>  order  Rungc  kutta  (Jill,  the  differential  equations  of  motion  are 
numerically  integrated. 


(4)  ( ONY I K I 1 or  a specified  sensor  location  ( latitude,  longitude  and  height)  the 
vehicle's  metric  data  and  their  time  rates  where  calculated  for  that  portion  of  its  flight  within  the 
sensor’s  2 0°  right  circular  cone 

( 5 ) PRIN 1 I he  time  histories  of  the  trajectories  are  presented  in  geocentric  and  topo- 
centric  coordinate  systems. 

((>>  BALCOI  Provides  a variable  ballistic  coefficient  generated  from  metric  data  for  each 
specific  vehicle.  I he  altitude  dependent  coefficients  were  extensively  employed  for  the  4°  of 
freedom  simulation  study. 


H )rth  Celestial  Pole 


Coordinate  System 

Figure  2-2  — Vehicle’s  Inertial  & Tropocentric  Coordinates 

2.1  Transfer  Maneuvers 


The  analysis  of  an  interception  with  a vehicle  requires  that  the  velocity  and  radius  vectors  of 
the  target  vehicle  be  matched.  The  analytical  analysis  was  reduced  to  the  problem  of  coplanar  orbital 
transfer.  The  transfer  maneuver  required  a definition  of  the  launch  timing  for  the  two  general  co- 
planar methods.  These  are: 

i)  Direct  ascent  to  the  synchronous  altitude  of  the  target. 

ii)  Launch  utilizing  parking  orbits. 

These  methods  neglect  perturbative  accelerations  and  assume  the  burning  time  of  the  inter- 
ceptor stages  to  be  infinitesimal. 
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For  the  computational  analysis  the  transfer  maneuver  is  assumed  to  be  elliptical  (coplanar 
and  non-coplanar)  and  the  target’s  orbit  to  be  circular.  Injection  is  assumed  to  occur  at  a point 
corresponding  to  the  intersection  of  the  transfer  elliptical  orbit  with  that  of  the  synchronous  target 
orbit.  The  relative  position  of  this  point  is  defined  with  respect  to  the  target’s  position  (at  the  time 
of  initiating  the  transfer  maneuver)  by  the  dwell  angle,  defined  in  Section  3.  A two-dimensional 
representation  of  the  above  described  scenarios;  methods  i and  ii  are  illustrated  in  Figures  2-3,  2-4 
and  2-5. 


Figure  2-3  - Geometry  of  the  Interceptor’s  Direct  Ascent  into  an  Elliptical 
Transfer  Maneuver 


I igure  2-4  Geometry  of  the  Interceptor's  Transfer  Maneuver  to  the  Higher  Altitude 
Synchronous  Satellite  Orbit 

2.2  Sensor  Locations 

The  previously  described  scenarios  can  be  exercised  against  a pre- 
selected set  of  sensor  locations. 


Figure  2-5  — Geometry  of  the  Interceptor's  Transfer  Maneuver  to  the  Lower 
Altitude  Synchronous  Satellite  Orbit 


SECTION  3 

DWELL  ANGLL/TKANSEEK  TIME 


\n  interceptor  in  launched  from  u specific  launch  site  on  a ballistic  trajectory  for  injecting 
into  a low-alliludc-dwcll  parking  orbit  at  a specific  inclination  A “time"  is  calculated  to  initiate  a 
maneuver  to  enable  the  interceptor  to  execute  a llobmann  elliptical  transfer  maneuver,  l or  the 
llohmann  transfer,  the  maneuver  is  initiated  at  the  mode  opposite  that  at  which  rendezvou  ■ is  to 
take  place.  I Inis  the  liinm  • must  be  precise,  that  is.  at  the  initiation  of  the  transfer  maneuver,  the 
s\  u hronous  ecpi.iii  i.il  N.iteiui  • target  mu  t he  at  a prescribed  orbital  position. 

I he  "time  at  which  this  translet  is  I "Mi.ited  is  defined  in  terms  of  an  angular  displace- 
ment. (-).  I his  angle,  defined  at  the  point  ot  departure  ol  the  interceptor,  is  a function  ol  the  ratio 
•l  ’.lie  orbital  radii  between  the  interceptor  and  the  target.  Defining  the  radial  distance  from  the 
earth's  center  to  the  interceptor’s  circular  orbit  and  the  distance  from  the  earth’s  center  to  the 
target  s circular  orbit  as  rj  and  rss.  respectively,  the  following  type  of  maneuver  is  presented. 

-'  l l Interceptor  Above  large!  I lie  interceptor  is  initially  on  a circular  orbit  with  a radius 
rj  greater  than  that  ot  the  target's  radius  rss.  i.e.,  the  interceptor  leads  the  target  by  the  angle  0 
(see  figure  3-1 ).  1 et  the  llight  time  ol  the  interceptor  s elliptical  transfer  maneuver  (apogee-perigee, 
7r>  initiated  by  an  impulsive  velocity  increment  at  apogee,  be  equivalent  to  the  target's  (light  time 
during  an  angular  displacement  (tt  + 0)  Since  the  target  exercises  uniform  angular  velocity,  the 
follow  ing  relationship  is  formed 


j * fM  ] 

l 'e‘  1 

1 I1  1 

where  I ^ is  the  period  of  the  i target’s)  circular  orbit  and  2 I e.‘  equal  to  the  period  of  the  elliptical 
orbit,  therefore. 


and  the  llight  time  ol  the  interceptor  is  reduced  (in  terms  of  the  variables  (ra.  rss»to  the  form: 


: r ' =•  a 


i..1 


i _ Jr 


ra  + rss 


1 1 ' - 


3/2 


(3) 


Substituting  equations  ( 2 1 and  1 3 ) into  equation  t I ) yields  the  following  expression  for  the  lead 
angle,  i.e..  the  interceptor  is  ‘ahead’  of  the  i iiget: 

<->  = * {,  l/2(ru,  rss  + l)|3/2-l| 

,or  ra  rss  I (4) 


I he  above  function  in  illustrated  in  figure  3 I 


INTERCEPTOR  LEADS  THE  TARGET 


9 = it  {[1/2(1  + ra/rgg)]3/2  - 1}  FOR  ra/rgg  * 1 

Figure  3-1  - Geometry  of  the  Relative  Positions  of  the  Target  and  an  Interceptor 
(Initially  Above  the  Target)  at  Initiation  of  the  External  Maneuver 

3- 1 -2.  Interceptor  Below  Target.  The  interceptor  is  initially  on  a circular  parking  orbit  with 
a radius  rp  less  than  that  of  the  synchronous  equatorial  satellite  radius  rss,  i.e.,  the  interceptor 
lags  the  target  by  the  angle  ©.  As  before,  the  “time”  to  execute  the  maneuver  is  given  by  this 
angle.  Employing  the  same  notations,  the  lag  angle  is  defined  in  the  following  manner.  Referring 
to  Figure  3-2,  the  following  relationship  is  formed: 


it  - 6 I _ I 2ff 

Te'  if  Tss 

where  Te‘  represents  the  time  the  target  uniformly  orbits  a (it  - 0)  angular  displacement. 


(5) 


3-2 


INTERCEPTOR  LAGS  THE  TARGET 


e =n{  1 - [1/2(1  + r /r„J]5/?}  FCR  r/r  < 1 

Uoo  P o o 

Figure  3-2  — Geometry  of  the  Relative  Positions  of  the  Target  and  an  Interceptor 
(Initially  Below  the  Target)  at  Initiation  of  the  Internal  Maneuver 

From  the  definitions  of  Tss  and  Te*  (equations  2,  3),  equation  (5)  may  he  written  in  terms  of 
the  variables  rp.  rss: 

0 = tr  {l  - [ l/2(rp/rss  + 1 ) 1 3/2  J 

for  Tp/rss  <1  (6) 

w here  the  position  of  the  interceptor  at  its  point  of  departure  is  equivalent  to  the  perigee  of  the 
elliptical  transfer  maneuver  ( rj  = rp)  and  its  point  of  rendezvous  equivalent  to  the  radius  of  the 
satellite/ target's  non-coplanar  synchronous  parking  orbit  (ra  = rss).  The  general  expession  for  the 
angular  displacement  of  an  interceptor  at  a radius  rj  is  given  by: 

0 = »r  {l  -M/2(ri/rss+  1)1  3/2 } (7) 

where  the  interceptor  may  lead  ( rj/ rss  > I ) or  lag  (rj/rss  < 1 ) the  synchronous  satellite/target.  The 
above  tuiu  tiwi. . graphically  illustrated  in  f igure  3-3  for  the  values  of  0 < rj/rS!S  < 3.0. 
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RATIO  OP  THE  INTERCEPTOR'S  RADIUS  TO  THE  SATELLITE  ORBIT  (R./R  ) 
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figure  I>well  Angle  vs.  Ratio  of  the  Interceptor's  Radius  to  the  Satellite  Orbit 


3.2.  DWFLL  ANGLF 


The  time  at  which  a transfer  maneuver  is  to  be  initiated  to  accomplish  a desired  operational 
mission  can  be  defined  in  terms  of  the  angular  separation  of  the  relative  geocentric  positions  of  the 
synchronous  satellite  (target)  and  the  interceptor.  I his  so-called  “dwell  angle  is  defined  lor  the 
final  elliptical  transfer  maneuver,  i.e.,  the  designed  point  ot  departure  from  the  interceptor’s 
parking  orbit  to  the  point  of  rendezvous  with  the  target’s  equatorial  parking  (synchronous)  orbit. 
These  parking  orbits  may  or  may  not  be  coplanar.  i he  dwell  angles  are  given  by  equations  (4)  and 
(6)  for  the  High-Altitude-Dwell  Scenario  and  the  Low-and  Mcdium-Altitude-Dwell  Scenarios, 
respectively.  A numerical  evaluation  of  these  equations  is  presented.  For  these  specified 
scenarios,  the  subsatellite  position  is  located  at  0°  latitude  and  80°  long- 
itude at  the  initialization  of  the  interceptor's  final  point  of  departure. 

3.2. 1 . Migh-Altitude-Dwcll.  This  particular  maneuver  defines  the  final  elliptical  transfer  stage 
of  the  bielliptic  transfer  scenario.  These  transfer  Hohmman  ellipses  (intermediate  and  final)  are 
unique  to  the  High-Altitude-Dwell  Scenario  and  are  illustrated  in  Figure  3-8.  Thus,  the  dwell  angle, 
given  by  equation  (4),  will  define  the  relative  geocentric  angular  separation  of  the  interceptor’s 
point  of  departure  (from  its  final  circular  parking  orbit)  and  the  synchronous  satellite  orbital 
position.  At  the  point  of  rendezvous  (t  = Te* ),  the  satellite  and  interceptor’s  total  angular  displace- 
ments are  (jt  + 0)  and  7r,  respectively.  Rewriting  equation  (4),  the  dwell  angle  for  the  interceptor’s 
final  elliptical  transfer  maneuver  becomes; 

0 = ir  {[l/2(raf/rss+  l)l3/2-l|  (7) 

where  raf  is  the  interceptor's  radius  at  the  apogee  of  the  final  elliptical  orbit  and  is  equivalent  to 
the  radius  of  the  interceptor’s  final  parking  orbit.  rcf.  From  Table  5-1.  the  nondimensionaliz.ed 
ratio; 

i _ i - (3*  1 93 10.54 1 1 3/6080.28  + 20925640,) 

raf' rss  rcf/rpf  13833730.1 

= 2.697461494  > I 

resulting  in  the  following  value  for  the  dwell  angle, 

0 = 272.469278  deg. 

The  functional  dependency  of  the  dwell  angle  on  the  non-dimensionalized  ratio  of  the  interceptor’s 
radius  at  the  point  of  departure  to  the  radius  of  the  satellite’s  synchronous  parking  orbit  is  illustrated 
in  Figure  3-3.  In  the  interest  of  mathematical  expediency,  the  interceptor’s  final  point  of  departure 
for  all  scenarios  occurred  at  the  intersection  of  the  orbital  planes  (i.e.,  at  the  ascending  node).  Thus, 
the  angular  separation  (dwell  angle)  was  easily  translated  into  a longitudinal  displacement.  The 
above  numerical  value  for0,  may  be  expressed  as  7. “5370722  W longitude  (interceptor’s  point  ol 
departure). 
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3 Low/Medium-Altitude-Dwell.  In  the  Low  and  Medium-Altitude-Dwell,  the  interceptor’s 

elliptical  transfer  maneuver  is  executed  from  a circular  parking  orbit  of  altitude  100  NM  and  3000 
NM.  respectively.  In  both  scenarios,  the  interceptor  (which  is  below  the  target)  lags  the  target,  i.e., 
the  interceptor  s and  target’s  total  angular  displacements  are  n and  (n  -0),  respectively  (see  f igure 
3-2).  Since  equation  ((>)  is  applicable  for  these  specific  scenarios,  the  numerical  value  of  0 follows. 

I or  the  Low-Altitude-Dwell  scenario.  Table  5-1  provides  the  values  of 

n/r  c = r /r  =r  /r  = 2.1533640  * I07  ft 

i ss  c ss  pi  a 1 3.8333730 j • io7  ft  01556644 

and.  for  the  Medium-Altitude-Dwelt  scenario, 

r./r  = r /,  = 3.9165640  ft  _n 

ss  J SS  13.83337301  ft  ~ 0-288  243  <1 

In  the  evaluation  of  equation  (6)  for  these  two  cases; 

0 =7r{l-|l/2(rj/rss+ l)|  3/2 } (8) 

0 = 100.936656  deg.  for  the  Low-Altitude-Dwell  and  0 = 87.502370  deg.  for  the  Medium-Altitude- 
Dwell.  With  the  satellite  located  at  its  previously  defined  longitudinal  position  (80,°0  H),  the 
angular  displacement  ot  the  interceptor  at  its  point  of  departure  is  given  by  20°. 936656  W 
longitude  and  7°. 502370  W longitude  for  the  Low-and  Medium-Altitude  dwell  scenarios,  respectively. 

3.3.  I LIGHT  T1MF 

With  the  definition  ol  the  target’s  period  in  a circular  orbit  about  a spherical  earth,  the  inter- 
ceptor s (light  time  during  the  elliptical  transfer  (Hohmann)  maneuver  may  be  rewritten  as: 


Te*  = ~7ihs'2  1 1/2 (rj/rss  + 1)1 


regardless  ot  whether  or  not  the  interceptor  leads  or  lags  the  target.  Figure  3-4  is  a plot  of  T'  as  a 
t unction  ot  the  parameter  ri/rss,  with  equation  (9)  applicable  in  this  case  for  the  entire  range  of 
ri/rss 

Another  factor  of  interest  is  the  period  of  the  elliptical  orbits  being  considered.  The  Keplcrian 
equation  tor  the  period  of  an  interceptor  in  an  elliptical  orbit  is  given  by: 

r - 2n  .,m 


where  the  semi-major  axis:  a = (ra  + rp)/2 


Tc  = l(rp  + ra)/2|: 

v in 
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Interceptor’s  Flight  Time  for  an  Elliptical  Transfer  Maneuver  vs.  Ratio 
of  the  Interceptor's  Radius  to  the  Satellite  Orbit 


which  may  be  reduced  in  terms  of  the  previously  employed  variables  (rj,  rss).  Therefore,  the  period 
of  the  interceptor  in  an  elliptical  orbit  is  given  by. 

Te 


Te  = 

where  rj/rss  < I should  the  interceptor  lag  the  satellite  target,  i.e.,  rj  = rp  and  rss  = ra,  and 
ri/rss  > I should  the  interceptor  lead  the  satellite  target,  i.e.,  rj  = ra  and  rss  = rp. 

With  the  definition  of  equation  (2),  the  above  equation  may  be  reduced  to  the  non- 
dimensionalized  orbital  period  of  the  interceptor: 

Te/Tss=[l/2(1 +ri/rss)|3/2  (14) 

Figure  3-5  is  a plot  of  the  above  equation  as  a function  of  the  parameters  r,/rss  for  the  entire  range 
of  interest. 


= “ l(ri  + rss)/2l 3/2  (12) 

P*  I* 

2v  rss3/2  11/2(1 +rj/rss)  1 3/2  (13) 


Another  expression  of  interest  may  be  obtained  by  rewriting  equation  ( 1 2)  as: 


Te  = A n3/:  [1/2(1  + rss/rj)l 3/2 
1 1/2 


(15) 


Since  the  penod  of  the  interceptor’s  circular  orbit  about  the  spherical  earth  may  be  written  as: 


7 = r 2'2 

C #i>/2  C 


(16) 


the  above  expression  may  be  reduced  to  the  nondimensionalized  form: 

Te/Tc  = [1/2(1+  rss/rj)l 3/2  (17) 

where  rj  = rc  and  rc  = ra  or  rc  = rp  if  the  interceptor  leads  or  lags  the  satellite/target,  respectively. 

For  this  particular  exercise,  the  target  is  located  in  an  equatorial  parking  orbit  at  synchronous 
altitude.  The  type  of  attack  modes  to  be  simulated  required  the  interceptor  to  initiate  a: 

(i)  Direct  Ascent  Maneuver 

(ii)  Low-Altitude-Dweil  Scenario 

(iii)  Medium-Altitude-Dwell  Scenario 

(iv)  High-Altitude-Dwell  Scenario 

I or  each  of  the  above  modes  (illustrated  in  Figures  3-6,  3-7,  and  3-8),  the  initial  state  vector  of  the 
interceptor  and  target  required  the  values  of  the  dwell  angle  0 and  the  interceptor’s  (light  time. 
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figure  3-6  Cieomctry  of  the  Interceptor's  Direct  Ascent  into  an  blliptical  Transfer 
Maneuver 

3.3.1  High-Altitude-Dwell  I light  I line  Analysis 

lo  calculate  the  interceptor's  (light  time  for  the  elliptical  transfer  maneuvers,  reference  is 
made  to  equation  (17). 

Rewriting  this  equation  for  the  intermediate  elliptical  transfer  maneuver; 

Tei  = Tcil  1/2(1  +rai/rci)|3/2  (18) 
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Figure  3-7  Geometry  of  the  Interceptor's  Transfer  Maneuver  to  the  Higher  Altitude 
Synchronous  Satellite  Orbit 

where: 

rej  = |ieriod  of  the  Hohmann  ellipse 

Ici  = period  of  the  initial  circular  parking  orbit 

and 

(raj/rcj)  = (raj/rpjl  = nondimensionalized  ratio  of  apogee/ perigee  parameter 

lo  calculate  the  value  of  the  period  for  the  initial  circular  parking  orbit,  with  an  altitude  of 
I DO  NM.  reference  is  made  to  equation  ( Hit.  Its  numerical  value  is: 

T ■ = -i!-  r,.;3/2  = 88.19808282  min.  (19) 

cl  pi/2  u 

where; 

rcj  = 2.1533640  * I07  Ft 
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F igure  3-8  - Geometry  of  the  Interceptor’s  Transfer  Maneuvers  from  a Low-Altitude- 
Dwell  to  a High-Altitude-Dwell  Parking  Orbit  to  a Final  Synchronous 
Altitude  Orbit 


Smte  the  apogee  radius  r;li  is  to  be  equivalent  to  three  times  the  altitude  of  the  synchronous  orbit 
its  numerical  value  is  determined  by 


and 

resulting  in 


rui  ~ 3 * 19310.541  13/6080.28  + 20925640.)  = 37.31499103  * I07  Ft. 
rai/rei  = 17.32869624 
Tei  = 88.19808282  (27.74293305)  min. 


2 


r 


or 

Tej  = 2446.873507  min. 


= period  of  the  intermediate  ellipse 

Tc'j  = Tei/ 2 = 1223.436754  min. 

= flight  time  of  the  intermediate  elliptical  transfer  maneuver 
( perigee  - to  - apogee ).  ( 20 ) 


The  elliptical  transfer  maneuver  pattern  is  illustrated  in  the  following  diagram. 


TT**F 

OF 

nooriFF 


The  interceptor’s  second  (or  final)  elliptical  transfer  maneuver  flight  time  can  be  determined  from 
equation  (II).  i.e. 


Tc  = 2 »r/p1/2  |(rp  + ra)/2| 3/2  (II) 

which  reduces  to  (for  the  specific  maneuver) 

Tf  = Tss  (1/2(1  +raf/rss)|3/2  (21) 

where 

I f * period  of  the  Hohinann  ellipse 

Tss  = period  of  the  synchronous  satellite  orbit 
and 

( ra|/rss>  = <raf/rpf*  = nondimensionali/ed  ratio  of  apogee/ perigee  parameters 


The  numerical  value  of  the  perigee,  equivalent  to  the  synchronous  orbit  radius,  was  an  additional 
constraint  imposed  Since  raj  = rcf  = raf  and  solving  equation  (2)  for  rss,  the  previously  defined 
nondimensionali/ed  ratio  takes  on  a value  of 


, _ 37.31499103  * I07  I t 

al  ss_  1 3.83337301  * I01  Ft 


(22) 


2.(»974(»I494 


resulting  in 


Tef  = (24  * 60)  (2.5 1 3682932)  min. 

Tef=  3619.703422  min. 

= period  of  the  final  elliptical  orbit 
and 

Tef  = Tef 12  = 1 809.85 1711  min.  ( 23 ) 

= fiight  time  of  the  final  elliptical 
transfer  maneuver  (apogee-to-perigec) 

rite  interceptor's  final  elliptical  night  maneuver  is  illustrated  in  the  following  diagram 


3.3.2.  Low/Medium- Altitude- Dwell  Flight  Time  Analysis 

Having  previously  defined  the  night  sequence  associated  with  the  Low-Altitude  and  Medium 
Altitude  scenarios,  the  interceptor’s  time  of  night  in  its  elliptical  transfer  maneuver  can  be 
numerically  evaluated  by  equation  (13).  This  equation  is  given  by: 


where 


rss/2  I l/2(r;/rss  + 1 )| 3/2 


Tc 


period  of  the  ellipse 
2ir  3/2 
Ml/J  C 


Tc  = period  of  the  interceptor’s  parking  orbit 


(13) 


(24) 
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and 


(25) 


Tss=2*rss3'V/2 

Tss  = period  of  the  synchronous  satellite’s 
equatorial  parking  orbit 


For  the  l.ow-Altitude-Dwell  scenario 


and 


and 


ri/rss  = 2. 1 533640  * I07  Ft/13.83337301  * I07  Ft 
= 0.155664421 

Tc  = 88.19808282  min.  parking  orbit 
Tss  = (24.  * 60.)  min. 

Tss  = 1440.0  min.  parking  orbit 


rhe  numerical  value  lor  the  period  of  the  ellipse  is  given  by 

Te=  1440.011/2  (1. 1 5566442 1)1 3/2  min. 
Te  = 632.5066752  min. 


The  flight  time  of  the  interceptor  from  its  point  of  departure  (t  = 0)  to  its  point  of  rendezvous 
(t  = Te* ) is  given  by 

Te1  = Te/2  = 316.2533376  min.  (26) 

as  verified  by  Figure  3-4. 


For  the  Medium-Altitude-Dwell  scenario 


rj/rss  = 165640  * 107  I t _ 0.283124296 

1 ss  13.83337301  * 107  Ft 

with  the  period  of  the  interceptor’s  parking  orbit  given  by 

Tc  = 216.341987  min. 

and  the  period  for  the  elliptical  transfer  orbit  by 

Te  = 1440.0  [ 1/2(1. 283 1 24296) | 3/2  min. 

Te  = 739.9810384  min. 

I hus  the  interceptor’s  (light  time  (from  perigee  to  apogee)  is  given  by 

Tt!  =Te/2  = 369.99052  min. 
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(27) 


Since  the  satellite  exhibits  a uniform  angular  motion  (d0/dt)  equivalent  to  the  earth’s  rotational 
velocity  (coe),  for  an  angular  displacement  of  (rr  - 0),  the  corresponding  flight  time  Tg , can  be 
simply  calculated  by  letting 


Te‘ 

Te 

for  the  Low-Altitude-ywell  and 

T* 

n 

for  the  Medium-Altitude-Dwell. 

For  the  High-Altitude-Dwell,  the  satellite’s  total  angular  displacement  is  given  by  (jt  + 0) 
resulting  in  a corresponding  flight  time  of; 

ti  _ ( 180.  + 272.4629278)  Deg 
e 360./ 1440  Peg/ Min 

T<!  = 1809.85171 1 Min 

as  verified  by  equation  (23). 

NOMENCLATURE 
Tc  = period  of  the  interceptor’s  circular  parking  orbit 
Tcj  = period  of  the  interceptor’s  initial  circular  parking  orbit 
Tcf  = period  of  the  interceptor’s  final  circular  parking  orbit 
Tsj  = period  of  the  interceptor’s  synchronous  parking  orbit 
T*  = period  of  the  equatorial  synchronous  satellite/target  parking  orbit 
Tg  = interceptor’s  maneuver  (elliptical  transfer)  flight  time 
I t!j  = interceptor’s  intermediate  maneuver  (elliptical  transfer)  flight  time 
I gt  = interceptor's  final  maneuver  (elliptical  transfer)  flight  time 
Te  = period  of  the  interceptor’s  elliptical  orbit 


_ Or  - 0)  _ Or  - 0) 


d0/dt 


co,. 


(28) 


= (180-  i UP, 936656)  Deg. 

360./ 1440.  Deg/M  in 

= 316.253376  Min. 


= (180. - 87.502370)  Deg 

360./ 1440  Deg/Min 

= 369.99052  Min. 
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NOMENCLATURE  (Continued) 

I ei  “ period  of  the  interceptor’s  intermediate  elliptical  orbit 
l et  = period  of  the  interceptor  final  elliptical  orbit 
rc  = radius  of  the  interceptor's  circular  parking  orbit 
rcj  = radius  of  the  interceptor’s  initial  circular  parking  orbit 
"r^f  = radius  of  the  interceptor’s  final  circular  parking  orbit 

rp  = radius  of  the  interceptor  at  the  perigee  of  the  elliptical  transfer  maneuver 

r pj  = radius  of  the  interceptor  at  the  perigee  ol  the  intermediate  elliptical  transfer  orbit 

raj  = radius  of  the  interceptor  at  the  apogee  of  the  intermediate  elliptical  transfer  orbit 

ra  = radius  of  the  interceptor  at  the  apogee  of  the  elliptical  transfer  maneuver 
rpf  = radius  of  the  interceptor  at  the  perigee  of  the  final  elliptical  transfer  maneuver 

raf  = radius  ot  the  interceptor  at  the  apogee  of  the  final  elliptical  transfer  maneuver 

rj  = interceptor’s  position  vector 

rss  = radius  of  the  satellitc/target’s  synchronous  parking  orbit 

V = interceptor’s  velocity  vector 

H = earth’s  gravitational  constant  = 

1.407639  * I016  M\ 
sec 

u>0  = earth's  rotational  velocity  = 

0.4375269048  * I O'2  rad/sec 
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SfcCTION  4 

DIRECT  ASCENT  ELLIPTIC  AL  TRANSFER  MANEUVER 
TO  SYNCHRONOUS  ORBIT 


I Methodology 

An  interceptor  is  launched  from  a specified  launch  site  on  a ballistic  trajectory  which 
intercepts  a satellite  in  a synchronous  equatorial  orbit.  The  direct  ascent  trajectories  are  defined 
at  the  termination  of  the  final-impulse  of  a multi-stage  booster  by  an  initial  state  vector  of  the 
interceptor.  By  the  definition  of  the  interceptor’s  initial  state  vector  (a  six-element  column 
matrix)  numerical  integration  of  the  differential  equations  of  motion  (previously  defined)  provides 
the  time  history  of  the  interceptor’s  elliptical  transfer  orbit  thus,  encompassing  the  free  flight 
maneuver  to  the  point  of  rendezvous,  as  illustrated  by  Figure  4-1.  In  addition,  the  Trajectory 
Simulation  Program  outputs  (when  appropriate)  a time-correlated  observation  data  set  (topocentric 
system)  for  a specified  sensor.  This  data  provides  an  observed  set  ol  metric  data  and  their  rates. 

This  data  is  pertinent  to  the  assignment  of  a figure-of-merit  to  each  of  the  specified  trajectory/ 
sensor  scenario. 

2 . Launch  Sites 

The  previously  described  can  be  exercised  against  a preselected  set  of 
sensor  locations. 
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Figure  4-1  Geometry  of  the  Interceptor’s  Direct  Ascent  Into  Elliptical  Transfer 
Maneuver 

3.  Initial  State  Vectors 


Hie  initial  state  vector  of  the  interceptor  and  the  satellite  target  are  defined  at  the  termination 
ol  the  final-impulse  (burnout).  Mathematically,  the  state  of  the  interceptor,  at  time  ti.  may  be 
represented  by  the  partitioned  matrix. 


I he  six-element  column  matrix  (previously  mentioned)  defines  the  interceptor’s  position 
coordinates  (latitude:  <J>.  longitude;  X,  altitude;  II)  and  velocity  coordinates  with  their  numerical 
values  (for  each  of  the  three  launch  sites)  specified  in  Tables  4- 1 and  4-2 . For  each  launch  site, 
the  components  of  the  interceptor’s  velocity  vector  </v7,  7,  (3)  are  defined  within  the  frame  of 
reference,  illustrated  in  Figure  4-3.  Since  the  timing  must  be  precise,  the  initial  state  vector  of  the 
interceptor  must  initiate  a collision-type  course  when  the  satellite  target  (designated  -SI ) is  at  a 
prescribed  position.  I he  initial  state  vector  of  SI  is  implicit  in  the  definition  of  a synchronous 
equatorial  satellite  at  a specific  longitude,  i.c., 

7 = 0 = w/2;  (yf  = (R,  +HlW,  = 10087.4548  FT/S 

II  = 19,  310.54  NM;  0 = 0.0°.  X = 80.0°l 


VERTICAL 


REENTRY  ANGLE  Y 
HEADING  ANGLE  B 
VELOCITY  MAGNITUDE  I VI 


Figure  4-  2 Interceptor’s  Intrinsic  Coordinate  System 
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4.  Sensor  Location/Ground  Track 


The  direct  ascent  scenario  can  be  exercised  against  a set  of  sensor 
locations  preselected  within  some  geographical  latitude  perimeter.  The 
geographical  locations  of  these  sensors  and  their  physical  locations  can 
be  illustrated  by  visibility  contours  plotted  upon  a Mercator  Projection 
Map.  These  contours  were  generated  by  programming  the  mathematical  equa- 
tions on  the  RADC/HP  9820A  Calculator  and  plotting  their  solutions  (by 
an  appropriate  coordinate  transformation)  with  the  RADC/HP9862A  Plotter. 

A description  of  these  equations  and  the  procedure  employed  are  described 
in  Reference  5.  It  must  be  emphasized  that  these  contours  are  restricted, 
but  by  no  means  limited,  to  a specific  set  of  parameters,  such  as  Rs 
(slant  range)  = 2000  NM  and  E (elevation  angle)  = 2.0  degrees. 

To  complement  the  sensor(s)  contour(s),  the  interceptor's  earth  trace 
(geographical  latitude  and  longitude)  from  the  three  (previously  discussed) 
launch  locations  can  be  plotted  by  the  9862A  Plotter  on  the  Mercator 
Projection  Map.  The  ground  track  data  are  an  output  of  the  Trajectory 
Simulation  Program. 
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SECTION  5 

LOW/MEI)IUM-ALTITUDE-l)WELL  ELLIPTIC  TRANSFER  MANEUVER 
TO  SYNCHRONOUS  ORBIT 

I . Methodology 

The  basic  dynamical  equations  describing  the  (light  sequence  of  the  interceptor  launched 
trom  a near  earth  circular  parking  orbit  dwell  into  a Hohinann  transfer  ellipse  for  rendezvous  with 
the  synchronous  orbit  satellite  target  can  be  defined  by  the  following  equations. 

i.  On  the  interceptor’s  inertial  velocity  in  a circular  parking  orbit  of  radius  "?c,  its  speed  can 
be  given  by 

ivci  = Mm  < n 

where  |VC|  = 25,567.4209  Ft/s 

•'or  |Tt-  | = 2.1533640*  107  Ft  (100  NM  altitude) 

and  |VC|  = 18,958.02238  Ft/s 

tor  |T^.  | = 3.9165640  * 107  Ft  (3000  NM  Altitude) 

with  1 | < |TSS|  ; i.c.,~r^s  is  equal  to  the  radius  of  the  synchronous  satellite  orbit. 

ii.  The  poirU  of  departure  of  the  interceptor  from  this  dwell  orbit  is  accomplished  by  an 
impulse  velocity  AV , initiating  a Hohmann  elliptical  transfer  maneuver  (refer  to  Figure  5-2). 

From  the  inverse-square  law  equation  (unit  mass) 


F = -pr/ 1 r I 3 = - pT/  l“r  I 


where  a = V d V / dr 


The  energy  integral  becomes: 


d V = - p 


resulting  in 


|VP|  2 =|va|  2 +2p  {l/|fp|  - l/|?a|} 

where  Vp,  Va  = interceptor’s  velocity  at  perigee  (t.,)  and  apogee  ( r^)  respectively.  Since  the 
mipulse  initiates  a point  of  departure  at  "rp  = ~fc  and  a point  of  revdezvous  at  "ra  ='rss,  where 


rc  < rss 


| Vp|  2 ~ | Va | 2 + -P  i l/l'rcl  “ l/|~?ss|} 
or  with  equation  ( I ),  the  above  reduces  to  the  form 

l^p|a-N,  + 2|v6|Mi-N/|tJ} 
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Since  Vp  > an  impulsive  velocity  is  to  be  added  at  "?p  to  initiate  the  transfer.  This  quantity 
can  be  determined  by  the  following  method.  From  the  law  of  the  conservation  of  angular  momentum 


r x V = constant 

(6) 

applied  at  rp  and  ra  (scalar  form) 

lral 

1 iVa  1=  irp  1 fvp  i=>i  va  1=  \ iTp  1 / iTa  1 f iVp  1 

(7) 

Equation  (5)  reduces  to: 

ivpi’  = 

firpi/iraif j ivp  12  +2  ivcij  |i  -iTci/irssif 

(8) 

where 

llVpl 

/ ivc  1 f = f i/2d+  nrc  1/  iTaDf  -'/2 

(9) 

and  _ 

ra  — 

rss 

Since  the  impulse  velocity  increment  required  to  initiate  this  specific  transfer  maneuver  is  given  by: 


IAV,  1=  IVPI-  lvcl=  {(  I Vp  I / IVCI)  - 1}  IVCI 
and  substituting  equation  (9)  for  the  value  of  ( i^p  I / l^c  1 ) results  in 


IAV!  1 

1>U 

II 

( 1 ra  1 / 1 > ^ 2 { 1 / 2 < •r’a  1 /■  • + 1)1 

-1/2.1] 

(10) 

or 

IAV]  1 

II 

<i 

n 

( 1 ra  1 / frc  1 ) 1 1/2(  lral/  1 r“c  1 ) ( lTa 

1/  ircl+  l)j 

I"1'2-!] 

Rewriting  equation  (10),  the  nondimcnsionaliz.ed  ratio  of  the  interceptor’s  initial  velocity  impulse 
to  its  initial  parking  orbit  velocity  is  graphically  illustrated  (by  Figure  5-1 ) as  a function  of  its  apogee 
(of  the  elliptical  transfer  maneuver)  to  the  parking  orbit  radius.  The  numerical  evaluation  of  the 
above  equation,  required  a value  of  the  apogee  radius  ?a  to  be  equivalent  to  the  radius  of  the 
synchronous  satellite  orbit.  Hence, 


lral=  ITSS  I = 13.8333730 


and  jlr,!/  lrc  if  = 


13.8333730  * 10^  Ft 
2. 1 533640  * lO7  Ft 


6.4240755 


107  Ft 


for  100  NM  parking  orbit  altitude  resulting  in 

IAV  | 1 = 8067.16206  Ft/s 
and  ltfcl  < lVpl=  33634.58296  Ft/s 
for  the  Low-Altitude-Dwell  scenario. 

The  Medium-Altitude-Dwell  scenario,  results  in 


{ 1 ^a 1 / ^c 1 f = 


13.8333730  * 107  Ft 

3.9 1 65640  * I07  Ft 


3.5320176 
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anil  the  incremental  velocity  numerical  value  (equation  10)  becomes;  !AV  jl  = 4710.63729  I t/s 
anil  IVCI  < IV,,  i = 23, 668. 65967  Ft/s.  The  direction  of  the  applied  impulse  is  illustrated  in  the 
following  Figure  5-2. 


figure  5-2  Geometry  of  the  Interceptor's  Transfer  Maneuver  to  the  Synchronous  Satellite  Orbit 


lii.  On  the  elliptical  transfer  orbit,  the  magnitude  of  the  interceptor’s  free-tliglit  velocity  can 
be  obtained  from  the  energy  integral,  i.e., 

IV  I = I Vp  I + 2fi  {l/ Irl- 1/ f?pl}  (111 

where  r is  the  interceptor’s  elliptical  position  vector. 

With  the  results  of  equation  (9)  and  letting  'rp  = fc.;  the  magnitude  of  the  interceptor’s  velocity  at 
the  apogee  position  is  given  by: 

• Va  I = I Vc  I {l/2(  lral/  frcl)(  lTa  I / lr“cl+  Ilf"1  ~ (12) 

Using  the  previously  defined  numerical  values,  the  above  equation  reduces  to. 

lVa  1=  5235.70788  Ft/s 

lor  the  Low-Altitude-Dwell  scenario  (i.e..  100  NM  orbit)  and  for  the  Medium-Altitude-Dwell: 


iv.  The  impulsive  velocity  increment  to  circularize  the  interceptor  at  apogee  (point  of 
rendezvous,  as  illustrated  in  Figure  5-2)  is 

AV2  = Vss-Va  (13) 

where  the  direction  of  the  impulse  velocity,  determined  by 

tfss  1 = /W  |fssf  = 10,087.4548  Ft/s  (14) 

required  an  addition  to  the  interceptor’s  “deficient”  circular  speed  (i.e.,  I Va I < IVSSI.)  The 
magnitude  of  the  impulse  velocity  can  be  shown  to  be  functionally  dependent  upon  the  ratio  of 
the  parking  orbits  radii.  With  the  scalar  form  of  equation  (13)  and  equation  (12), 


IAV2I=IVSSI-I^al 

(13.1) 

and 

l 

IAV2I  = IVssI-IVcI  {l/2(  lral/  frcl)(  lral/  lTcl+  1)}  '1/2 

(15) 

Since 

IVCI  =/„/  l?ci(  l?ssl/  lfssl)  1/2  = IVSSI(  l?Ml/  r?el ) 1/2 

(16) 

and  fss 

= r^,  equation  (15)  reduces  to 

IAV2I=  IVSSI  {l-(l/2(  lrssl/lrcl+  l)]~l/2} 

(17) 

with  a numerical  value  of 


IAV2I  = 4851.74692  Ft/s 
for  the  Low-Altitudc-JJwell  scenario,  and 

IAV2I  = 3386.28219  Ft/s 

for  the  Medium-Altitude-Dwell. 

Rewriting  equation  (17)  in  a nondimensionalized  form,  i.e., 

(lAV2l/lVssl)={l-|l/2(lrssl/lTcl+  l))',/2f  (18) 

a plot  of  ( I AV  t I / 1 Vssi ) vs  the  ratio  of  the  circular  parking  orbits  radii  ( I rss  I / frc  I > 1 ) is 
illustrated  in  Figure  5-3. 

In  the  performance  of  these  specific  elliptical  transfer  maneuvers, 
the  "time"  associated  with  the  interceptor/satellite  (target)  relative 
positions  is  calculated  in  Section  3.  A graphical  display  of  the  inter- 
ceptor's velocity  and  altitude  at  these  specified  events  for  the  two 
scenarios  is  illustrated  in  Figure  5-4. 
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SECTION  6 

HIGH  ALTITUDE  DWELL  BIELLIPTIC  TRANSFER  MANEUVER 
TO  SYNCHRONOUS  ORBIT* 


I.  Methodology 

The  (light  sequence  of  the  interceptor  illustrated  in  Figure  6-1,  consists  of  a power  (light 
trom  a specific  launch  s^te  to  injection  into  a low-aititude-dwell  circular  parking  orbit  of  radius 
rcj  and  orbital  velocity  Vcj.  At  a specific  launch  “time”,  defined  in  terms  of  an  angular  displacement, 
an  applied  impulsive  velocity  increment,  defined  below  by  equation  (2),  injects  the  interceptor  in 
an  Hohmann  transfer  ellipse  of  an  apogee  radius  "faj,  illustrated  in  Figure  6-2.  The  point  of  depar- 
ture of  the  interceptor  from  its  parking  orbit  constitutes  the  perigee  radius  "fpj  ( = ?cj)  of  the 
elliptical  transfer  maneuver  of  the  intermediate  (first)  orbit.  Note  that  the  establishment  of  an 
elliptic  orbit  of  a given  parameter  was  by  means  of  an  impulsive  velocity  AV, . As  the  radius  of 
the  interceptor  approaches  the  value  of  the  apogee  radius,  a successive  stage  of  the  interceptor’s 
booster  is  ignited  to  provide  an  impulsive  incremental  velocity  AV2 . This  impulse  circularizes  the 
interceptor  orbit  at  apogee  (infinitesimal  burning  time),  with  a radius  ’rcf.  Since  the  velocity 
vector  of  the  interceptor  approaching  apogee  has  a “deficient”  circular  speed,  the  applied  velocity 
increment  provided  the  corrective  maneuver  i.e. , an  acceleration.  With  this  maneuver  the  inter- 
ceptor is  now  in  a High-Altitude-Dwell  circular  parking  orbit  of  radius  ~rcj\ 

Since  “time”,  at  which  these  transfers  are  to  be  initiated,  is  very  crucial,  the  final  elliptical 
transfer  maneuver  must  be  executed  at  the  node  opposite  the  point  of  rendezvous  (perigee),  to 
provide  the  necessary  correlation  with  the  target’s  (light  time  and  hence  its  angular  displacement 
with  respect  to  the  point  of  departure.  The  Hohmann  transfer  is  initiated  at  the  apogee  (point  of 
departure)  illustrated  in  Figure  6-4,  by  the  execution  of  an  impulsive  velocity  increment  AV,, 
in  the  direction  opposite  the  velocity  vector  since  the  interceptor  has  an  “excessive”  apogee 
velocity  - 

I he  equations  used  for  each  of  the  transfer  maneuvers  illustrated  in  Figures  6-2,  6-4  are: 

i.)  On  the  near-earth  (initial)  circular  parking  orbit  of  radius  Tcj,  the  speed  of  the  velocity 
vector  is  given  by 

IV  -1=  /jT  - A0675004  * I O'”'  FT3/MINJ  ... 

C1  /rci  "/ 27153364  MO7  FT  (l) 

lVcil=  25567.4209  FT/S 
for  a parking  orbit  altitude  of  100  NM. 


•Bielliptic  Transfer:  The  apogee  altitude  of  the  intermediate  (first)  transfer  ellipse  is  set  at  an 
altitude  three  times  the  altitude  ot  the  24-hour  altitude.  The  final  ellipse  is  then  required  to 
have  its  perigee  at  the  24-hour  orbit  altitude. 
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Figure  6-2  - Interceptor's  Intermediate  Elliptical  Transfer  Orbit  from  an  Initial-to- 
Final  (Desired)  Circular  Parking  Orbit 

ii.)  The  impulse  velocity  increment  required  to  initiate  the  intermediate  elliptical  orbit  is 
given  by; 

= ^pi  ” ^ci 


where 

IAV,  1=  lVcil[(rai/rci)‘'2  + l) 

} ']  • 

(2 

(at  the  perigee) 

For  the  numerical  evaluation  of  the  above  equation,  an  additional  constraint  required  the  value 
of  the  apogee  radius  - ?ai,  to  the  equivalent  to  three  times  the  altitude  of  the  equatorial  synchronous 
satellite  orbit.  Therefore; 

^ai1  = (3  * 19310.541 13/6080.28+  20925640.)=  373,149,910.3  FT 
and 

<lr.il/lfdl)-  1 7.3286964  2=<rcf/rpi)  = (raf/rpi) 
resulting  in  I AV,  I = 9590. 1 674  FT/S 

I he  direction  of  the  applied  impulse  is  illustrated  in  Figure  6-2.  Rewriting  equation  (2)  in  the  form 

o»  the  nondimensionalized  ratio  of  the  interceptor’s  velocity  impulse  to  its  initial  parking  orbit 
velocity,  the  following  equation  is  graphically  illustrated  in  Figure  6-3. 


iii.)  On  the  intermediate  ellipticai  transfer  orbit,  the  magnitude  of  the  interceptor’s  velocity, 
functionally  dependent  on  the  ratio  of  the  apogee  to  initial  circular  parking  orbit  radii,  is  given  by; 


(at  the  perigee) 

And  from  the  conservation  of  angular  momentum,  the  magnitude  of  the  apogee  velocity  may  be 
expressed  by; 

l^ai'  = ^piK^pil/lrai1)-  (4) 

Substitution  of  equation  (3)  for  I Vpj  I reduces  the  above  expression  to 

*-■- 'Mi  ($)(£♦■)}'“ 

(at  the  apogee) 


Figure  6-4  Interceptor’s  Final  Elliptical  Transfer  Orbit  to  the  Synchronous 
Satellite  Target  Orbit 


With  the  numerical  values,  previously  defined,  equations (3 1,  (5)  become 


I Vpi  I = 35  1 57.5883 1 IT/S 

ltfail=  2028.865  I 5 I T/S 

(in  a C C W direction,  as  illustrated  on 
page  3- 1 3 ). 


iv. ) Die  impulsive  velocity  increment  required  to  circularize  the  interceptor’s  final  parking 
orbit  is 


AV2=Vcf-Vai 

for  lVcfl  > IVall 

with 


lV,.fl  = = 6141.91744  FT/S 

/ rcf 

and  the  previously  calculated  value  for  lVajl , results  in 

lAVj  I = 41  13.05229  FT/S 


(<>) 
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applied  in  the  direction  of  motion  (clockwise  illustrated  in  Figure  6-2)  due  to  a “deficient” 
circular  speed,  exemplified  by  equation  (6). 

v.)  The  magnitude  of  the  impulsive  velocity  increment  required  to  initiate  the  final  elliptical 
transfer  orbit  can  be  determined  from  the  following  method.  The  speed  of  the  interceptor’s  velocity 
at  apogee  is  given  by 


(8) 


(at  the  apogee) 


where  the  numerical  ratio  of  the  interceptor’s  final  circular  parking  orbit-to-the  target’s  synchronous 
parking  orbit  radii  is  given  by 


resulting  in 


(a  c:  h 


373 1499 10. 3\ 
138333730.1/ 


= 2.697461494 


IVafl  = 4517.17736  FT/S 


Since  the  speed  of  the  interceptor’s  velocity  in  the  final  parking  orbit  lVcfl>  lVaj  l . the  “excessive” 
apogee  velocity  must  be  reduced  at  ?af  via  the  equation 

AV,  =Vcf-Vaf 
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where 


IAV,I  = lVcfl  - lVafl 

= (>141.91744  - 4517.17736  (9) 

lAVjl  = 1624.74008  FT/S 
(counter  to  the  interceptor’s  direction). 

vi.)  The  magnitude  of  the  impulsive  velocity  increment  required  to  circularize  the  inter- 
ceptor into  the  target’s  equatorial  synchronous  parking  orbit  is 

iav4i  = ivpfi  - lvssl  (10) 

or 

(at  the  perigee) 

for  rcf  > tss 

From  the  previously  defined  numerical  values  of  I Vssl  and  the  radii  ratio;  (rcf/r^s),  the  interceptor’s 
“excessive”  circular  velocity  requirement  is  reduced  by 

IAVJ  = 12184.91  199- 10087.4548  FT/S 
IAVJ  = 2097.45719  FT/S 

Since  the  flight  of  the  interceptor  on  a Hohmann  elliptical  transfer  maneuver  is  on  a collision-type 
course,  it  need  not  circularize  at  the  synchronous  satellite  target  orbit.  This  concept  reduces  the 
total  required  incremental  velocities  (i.e.  AV4  -*  0)  for  the  sequence  of  maneuvers. 

Rewriting  equation  ( 1 I ),  a plot  of  the  interceptor’s  nondimensionalized  final  impulse  required 
to  circularize  its  perigee  velocity.  at  synchronous  altitude  is  given  by  Figure  6-5. 

In  summary,  the  interceptor’s  altitude  and  inertial  velocity  for  each  of  its  specified  events 
during  the  transfer  maneuvers  are  presented  in  Figure  6-6. 


With  the  above-defined  numerical  values,  the  simulation  of  the  interceptor’s  final  elliptical 
transfer  maneuver  to  the  point  of  rendezvous  (with  the  target  S 1 ) was  accomplished  by  the 
previously  defined  Trajectory  Simulation  Program.  This  simulation  (of  the  interceptor’s  tlight 
time  of  ~ 1800  min)  was  restricted  to  the  single,  zero-inclination  elliptical  transfer  orbit,  and 
employed  the  data  (flight  time,  dwell  angle)  derived  in  Section  3. 
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INTERCEPTOR’S  NONE MENS ION ALI ZED  FINAL  IMPULSE  TO  CIRCULARIZE  AT  SYNCHRONOUS  ALTITUDE  ( aV 


Figure  <>-S  K;itio  of  Interceptor's  Final  Parking-to-Synchronous  Orbits  Radii 
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(Applicable  to  the  liigh-Altitude-Dwell  Bielliptic  I ransfer  Maneuver) 

Vpj  = velocity  of  the  interceptor  ai  the  perigee  of  the  intermediate  elliptical  transfer  orbit. 

\7al  = velocity  of  the  interceptor  at  the  apogee  of  the  intermediate  elliptical  transfer  orbit. 

Ccl  - velocity  of  the  interceptor  in  the  initial  circular  parking  orbit. 

Vcf  = velocity  of  the  interceptor  in  the  final  circular  parking  orbit. 

— » 

Vss  = velocity  of  the  target  in  the  equatorial  synchronous  parking  orbit. 

Vpj'  = velocity  of  the  interceptor  at  the  perigee  of  the  final  elliptical  transfer  orbit. 

-Jk 

Va|  = velocity  of  the  interceptor  at  the  apogee  of  the  final  elliptical  transfer  orbit. 
fcj  = radius  of  the  interceptor’s  initial  circular  parking  orbit. 


■?cf  = radius  of  the  interceptor’s  final  circular  parking  orbit. 
r“sS  = radius  of  the  target’s  synchronous  parking  orbit. 


AV, 
r 

AVj 


AV3 

AV4 


= Cp!  - Vcj  = difference  in  the  interceptor’s  perigee  velocity  and  its  circular  velocity  for 
the  intermediate  (initial)  transfer  orbit. 

= Vcf  - Vai  = difference  in  the  interceptor’s  velocity  in  the  final  circular  orbit  and  the 
apogee  velocity  in  the  initial  transfer  orbit 

= Vcf  - Vaf  = difference  in  the  interceptor’s  velocity  in  the  final  circular  orbit  and  the 
apogee  velocity  in  the  final  transfer  orbit 

fe  k 

= Vpf-  Vjj  = difference  in  the  interceptor’s  perigee  velocity  in  the  final  transfer  orbit 
and  the  circular  velocity  in  the  desired  (synchronous)  orbit. 


i 
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APPENDIX  A 

INERTIAL  GEOCENTRIC  ANGULAR  SEPARATION  OF  THE  DIRECT  ASCENT 
INTERCEPTOR/TARGET  SCENARIO 

For  the  purpose  of  early  warning,  a time  related  angular  separation  of  specific  interceptor 
launch/target  location  scenarios  were  calculated  by  the  Trajectory  Simulation  Program  and  plotted 
with  the  HP  9862A  plotter.  These  plots,  the  longitudinal,  and  inertial  geocentric  angular  separation 
include  the  Direct  Ascent  interceptor  maneuver  to  the  equatorial  synchronous  satellite/target  and 
are  illustrated  in  Figures  A- 1 and  A-2  for  the  three  launch  locations.  These  angular  separations  are 
calculated  during  the  interceptor’s  total  flight  time  ( ss  149.4  minutes)  referenced  to  the  zero  time 
at  collision.  Specifying  the  interceptor  burnout  state  vector  and  the  target’s  time  corresponding 
state  vector,  the  time-dependent  angular  separation  was  calculated  by  the  Trajectory  Simulation 
Program  and  the  results  plotted  by  the  HP  9862A  plotter.  The  mathematical  derivation  of  the 
angular  displacement  is  given  by 

a(t)  = COS-1  | COS  (AX)  COS  0j  COS  0,  + SIN  0j  SIN  0,1  (A-l) 

where 

AX  = (Xj(t)  - Xt(t)|  (A-2) 

AX  = longitudinal  separation  of  the 
and  interceptor/ target 

0,<t)  = interceptor’s  geocentric  latitude 

0(<U=  target’s  geocentric  latitude 

For  the  equatorial  synchronous  satellite  target  0 = 0.0.  Deg  I and  the  above  expression  reduces  to. 

o<t>  = COS-'  | COS  (AX)  COS  0t)  > A-3) 

The  interceptor/ target  geometry  is  illustrated  by  the  following  figure. 

Reference  to  Figure  A-2  illustrates  that  the  inertial  geocentric  angular  separation  at  burnout 
time  is  functionally  dependent  on  the  location  of  the  specific  launch  site  and  may  vary  from 
approximately  45  degrees  to  70  degrees.  Initially,  the  angular  (closure)  rate  is  relatively  high  and 
approaches  a moderate  constant  value  approximately  100  minutes  prior  to  the  intercept. 

The  longitudinal  angular  displacement,  illustrated  by  Figure  A-l  is  initially  negative,  illustrating 
that  the  target’s  longitudinal  position  is  further  I ast  as  compared  to  the  interceptor  position.  The 
angular  rate  is  quite  high  and  does  not  approach  a constant  value  until  about  100  minutes  prior  to 
collision. 


A-l 


Figure  A-l.  Longitudinal  Separation  of  the  Interceptor/Target  Prior  to  Collision 
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Figure  A-2  Inertial  Geocentric  Angular  Separation  of  the  Interceptor  Target  Prior  to  Collision 


APPENDIX  li 

INTERCEPTOR  CO-PLANAR  CIRCULAR/ELL1P1 ICAL  ORBIT  ANGULAR  SEPARATION 


In  order  to  lx-  able  to  detect  the  time  of  departure  of  an  interceptor  from  its  circular  parking 
orbit  onto  a Hohmann  elliptical  transfer  maneuver  to  the  point  of  rendezvous,  a study  was  made  to 
determine  the  time-correlated  angular  separation. 

With  the  Trajectory  Simulation  Computer  Program,  a circular  equatorial  parking  orbit  and  a 
coplanar  elliptical  transfer  orbit  were  simulated  for  the  three  altitude-dwell  scenarios.  I he  initial 
position  and  velocity  of  the  interceptor  in  a circular  and  elliptical  motion  are  given  in  Table  B-l . 

For  each  simulated  scenario,  the  interceptor’s  time-dependent  longitudinal  position  (incre- 
mented in  0. 1 minute  intervals)  in  the  circular  parking  orbit  was  correlated  with  its  elliptical  transfer 
maneuver  ( time-sequenced  with-respect-to  the  point  of  final  departure,  illustrated  in  Figure  B-l ). 

The  elliptical  transfer  maneuver  for  the  Low/Medium-Altitude-Dwell  Scenarios  occurred  at  the 
perigee  and  the  interceptor’s  longitudinal  angular  displacement  is  equivalent  to  the  true  anomaly 
whereas  the  High-Altitude-Dwell  Scenario  required  the  interceptor's  elliptical  transfer  maneuver  to 
be  initiated  at  apogee,  (equivalent  to  a true  anomaly  value  of  180.  degrees). 


Figure  B-l  Interceptor’s  Angular  Displacement 

In  each  ol  the  Altitude-Dwell  Scenarios,  these  longitudinal  angular  separations  (AX)  were 
plotted  as  a function  of  the  transfer  initiation  time  (T).  In  Figure  B-2,  the  Low-Altitude-pwell  » 

Scenario  exhibits  a greater  time  rate-of-change  of  angular  separation  (AX/At)  with  a value  of  6.° 
m approximately  5.  minutes  whereas  the  Medium-Altitude-Dwell  Scenario  required  approximately 
21.  minutes.  In  Figure  B-3,  the  High-Altitude-Dwell  Scenario  with  the  least  rate-of-change  required 
approximj  telv  J5().  minutes  of  elliptical  flight  time  prior  to  achieving  this  value.  These  values 
represent  1.0%,  5.7%  and  24 .9%  of  the  total  elliptical  transfer  maneuver  (light  time  for  the  I ow. 
Medium  and  High-Altitude  Dwell  Scenarios  respectively. 


B-l 


TABLE  B-1 

INITIAL  STATE  VECTOR  OF  THE  INTERCEPTOR  AT  TRANSFER  TIME 


MEDIUM  ALTITUDE  DWELL 
(3000  NM) 


TIME  FROM  INITIATION  OF  FINAL  ELLIPTICAL  TRANSFER  MANEUVER  T(MIN) 


Figure  B- 1 High- Altitude-Dwell  Angular  Separation 


APPENDIX  C 

ORBITAL  TRANSFER  ENERGY  REQUIREMENT 

The  total  velocity  increment  A7r,  required  by  an  interceptor  to  execute  the  bielliptic 
transfer  to  a coplanar  synchronous  satellite  target  orbit,  as  demonstrated  by  the  High-Altitude- 
Dwell  Scenario  is  defined  by: 

i.  difference  in  the  initial  circular  parking  orbit  velocity  and  the  initial  perigee  velocity 
for  the  intermediate  (initial)  elliptical  transfer  maneuver  A-?, , plus 

ii.  difference  in  the  initial  apogee  velocity  for  the  initial  elliptical  orbit  anil  the  interceptor's 
circular  velocity  in  its  parking  orbit  (final)  A t>.  plus 

iii.  difference  in  the  circular  velocity  and  the  apogee  velocity  for  the  final  elliptical  transfer 
maneuver-Al?3  plus 

iv.  difference  in  the  perigee  velocity  in  the  final  elliptical  orbit  and  the  circular  velocity  of 
the  desired  orbit  Av4 ; i.e., 


Avy  =A’?1  +AT?>2  +AT?3  + A"?4  (C'-l) 

A four-velocity  increment  transfer  between  the  circular  (initial)  and  desired  circular  orbits. 

IAI?tI=  lAni  1+  lAi'jl  + lAi^l  + I A"?4 
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where  Y s ( rss/rcf) 


To  evaluate:  lAT^  I = ( - IT^SI  > 
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for  X > I , Y < 1 

Since  the  Ilohmann  transfer  ellipse  has  long  been  recognized  as  the  minimum  energy  for  a two- 
impulse  velocity  increment  between  coplanar  circular  orbits,  an  analysis  to  determine  whether  this 
is  the  most  efficient  means  of  transfer,  as  compared  to  the  four-velocity  increment  is  required.  I or 
the  two-velocity  increment,  the  transfer  between  the  initial  and  desired  circular  orbits,  equation 
reduces  to; 

=r  p*  -ii  + i?ssi  f,-/r  ((-25) 

Itv-ji  II  L/l  + x J 'K'ii  l / x /,  + x J 


a _ 


where  X = rcf/rd  = rss/rci>  I 
and  Ai^  I = IA1?4  i = 0 since  Y = I 

(liven  the  following  numerical  values  for  the  four-impulse  transfer,  equation  is  evaluated. 


X = (rcf/rci)=  373149910.3/21533640.0  (C-26) 

X = 1 7.32866942 

Y = CWrcf)  = 138333730. 1/373149910.3  (C-27) 

Y = 0.3707189156 

t^cf/^ci) = 6 14 1.91 744 1/25567.4209  (C-28) 

= 0.2402243646 


<"ss/*Vi>=  10087.4548/25567.4209 
= 0.3945433073 


resulting  in; 


(C-29) 


0.375074432  + 0.6696690938  + 
0.635472812  + 0.0820369968 


/IAP*t  IN 

I J-  0.681 5288426  total  energy  requirement  for  the  four-impulse  coplanar  transfer 
ci  ' between  circular  orbits  T\.j  and  (C-301 


For  the  two-impulse  energy  requirement  between  the  coplanar  initial  and  desired  circular  parking 
orbits,  equation  is  evaluated  with  the  following  numerical  data. 


Let  X =(rss/rcj)=  138333730.1/21533640. 


((  -3 1 ) 


X - (>.424075544 


Y = (rss/rss)=  1.0  (C-3  2 ) 

(•'ssM.il = 10087.4548/25567.4209  (C-33) 


= 0.3945433073 

resulting  in; 

/lAtf,  r , 

I , I 0.505287921  7 total  energy  required  tor  the  two-impulse  coplanar 
Mi1  ' transfer  maneuver. 


(C-3  4) 


f hus,  the  liohmann  two-impulse  transfer  is  the  most  economical  approach  for  this  specific  data  set. 
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APPENDIX  D 

AN  ANALYTICAL  SOLUTION  TO  THE  THREE 
DIMENSIONAL  REPRESENTATION  OF  THE  BALLISTIC  MISSILl 

The  motion  of  a ballistic  missile,  after  thrust  termination,  is  described  by  means  of  a non- 
rotating  spherical  polar  coordinate  system  (r,  0,  X)  with  origin  at  the  earth  s center  and  the  polar 
axis  coincident  with  that  of  the  earth. 


The  Lagrangian  in  this  system  is 


L = T V 


where  T = kinetic  energy 
T = m(v  ■ v)  / 2 

and  V=  V ( r)  = gm/r  (potential  energy) 

From  this,  the  equations  of  motion  are  derived  from  the  Euler-Lagrange  equation  in  general 
coordinates 


where  Oj  represents  the  nongravitational  forces. 

Since  the  velocity,  in  the  nonrotating  system,  may  be  expressed  in  terms  of  the  spherical 
coordinates  as 

v = r?+  r6  0 + (rCOS0)X  X 

and  vJ  = "v  ■ v = r2  + <rO)2+(rCOS  0)1  X2 

L = l r2  + r2  02  + (rCOS  0)2  X2  ] + 

Letting  q,  and  qj  be  equivalent  to  the  spherical  polar  coordinates  and  their  rates  respectively 
equation  D-3reduces  to  the  following  equations  of  motion  i.e., 

^7  - r 02-(rCOS2  0)X2  + ^ = Q,  \ 


( r2  0 ) + 


^7  SIN  0^  X2  =< 


l(r2COS2  0)  X|  =g., 


D-l 


The  initial  state  vector  of  the  missile  may  be  written  as 


dr 

— v 

dt 

O 

- v o 

d0 

= v° 

dt 

o 

” 'o 

dX 

= v« 

dt 

o 

fo 

COS  7 

SIN  7 COS  0 

SIN  7 SIN  p 
COS  0O 


Figure  D-l  - Spherical  Polar  Coordinates  of  the  Velocity  Vector 

I he  above  2nd  order  non-linear  differential  equations  of  motion  were  numerically  integrated  by 
the  4th  order  Runge  Kutta-Gill  Method  (1,5).  Since  the  Q’s  will  give  rise  to  a torque  and  an 
analytical  solution  is  required  to  establish  the  initial  state  vector  on  the  missile,  the  right  hand 
side  is  set  equal  to  zero.  This  results  in  the  angular  momentum  (time  rate  of  change  of  torque) 
being  conserved  and  the  Lagrangian  equation  can  be  rewritten  as 

(r2C()S 2 0)  X = constant 

and  the  remaining  equations  reduce  to 

d\r  -r02  - (rC'OS20 ) X2  + H =Q 
dt2  r2 

d ( r2  0 ) + r*  SIN  2 0)  X2  = O 
ut  2 


(D-7) 


(D-8) 

(D-9) 
(D- 10) 


D-2 


Since  the  angular  momentum  is  conserved,  the  expression,  as  given  by  equation  (D-8)  may  be 
rewritten  as 


H = (r2(/OS20)  dX/dt 
H = r0  vQ  COS  0o  SIN  7 SIN  0 

To  eliminate  the  time  dependency  of  the  Lagrangian  equations  tie.,  r(t),  0(1).  X<tt).  the  expression 
for  angular  momentum  provides  the  functional  dependency  of  the  differential  change  in  time,  dt 
to  the  corresponding  change  in  dX. 


CD-11) 


d - dX  d _ H _d_  fD-12) 

dt  " dt  dX  r2COSf0  dX 


Rewriting  equation  (D-10) 


Since 


d (r2  ())  + (r2SIN  OC'OSO)  X2  = O 
dt 


dX  _ H 
dt  r2  COS2  0 


d 

dt 


(r2  0)  + 


II2  tan  0 
r2(  OS20 


= 0 


CD-13) 


CD-14) 


CD- 15) 


I irsi  term  of  the  above  expression  (equation  (D- 15)  may  be  reduced  to 

({2  d 0\  It  d / II  m\ 

\ dt  J r2(  OS2«  d*  \tOS20  / 


d 

dt 


resulting  in 


since 


d 

/si  C20  dd 

dX 

y dX 

d 

tan  0 = SI  C' 

dX 

d 

( d a\ 

dX 

(^dX  Un0) 

) 


+ Ian  0=0 


d0 

dX 


dX2 


tan  0 


Solution  to  the  2nd  order  ordinary  differential  equation  (0-17) 

d2 


dX2 


tan  0 + tan  0 = 0 


Let  S = tan  0 and  ds/dX  = T 

Rewrite  the  above  equation  as  a first  order  equation. 


d2 S = d l (S)  „ 
d'X*  dX 


S)=  / d 1 \ / dS \ 

ydS ) l dX  J 


= TdT 
dS 


TdI  + s-o 

dS 

Solution  of  the  above  first  order  equation,  yields  a constant  of  integration  given  by: 


C,  =T2  +S2  = tan  0^  + tan2  0 

d tan  0 = SEC'2  0 (d0/dX) 
dA 


where 

a,'d  = (dO/dt)  / (dX/dt) 

I rom  the  initial  boundary  conditions  given  by  equation  (0-7) 


dX 


tan  0 


= (ct1&  \ 
o \<™o) 


and 


Cj  = (ctn2  P + SIN2  0O)  1/2 /COS  0o 


(0-16) 

(0-17) 

(0- 18) 
(D- 19) 

(D-20) 
(0-21 ) 

(0-22) 


(0-23) 

(0-24) 

(0-25) 

(0-26) 

(0-27) 


0-4 


Rewriting  equation  (D- 21) 


(0-28) 


and  integrating  again,  results  in 


tanfl  = ye?SIN(X  + C,) 
and 

SEC  0 = 11  +C,SIN2(X  + C, )j  1/2 


(0-29) 

(0-30) 


with  C, , C2  as  arbitrary  constants  of  integrations.  Differentiating  of  equation  (D-29)  and  imposing 
the  initial  set  of  boundary  values  on  the  two  equations,  results  in 


tan  0O  = ycTslNiXo  + C, ) 


(0-31) 


and 


d 

dX 


tan  0 


o 


yc7cos(X0+c,) 


(D- 32 ) 


Substituting  equation(D-26)  into  the  above  differential  equation  and  dividing  into  equation  (D-31)  , 
results  in  an  expression  representing  the  intersection  of  a trajectory  plane  with  a sphere  (see 
Figure  D-3. 

tan  (X0  + C, ) =SIN0o  tan  0 /-p 


EQUATOR 


POLAR  AXIS 


f igure  D-3  Unperturbed  Planar  Trajectory  on  a Spherical  Earth 


lo  impose  an  additional  constraint  upon  the  initial  state  vector,  let  X()  - ()  and  with  the  substitution 
ol  equation (D-26)in  equation  (D-32), the  resulting  equation  and  equation  (D-31)  reduces  to: 


w 

& 

ii 

so. 

(D-34) 

COS0o  * 

tan  0O  = yf^SlNC  , 

(D-3S) 

tan  0COS0o  =S1NX  ctn  P +COS  X SIN  0O 

tan  ( , = tan  OqUcXu  0/COS  0o) 

(D-3b) 

as  illustrated  by  Figure  l)-3. 


Also  verification  of  equation  (D-27))can  be  obtained  by  the  squaring  of  equations  (D-34) 
and  (D-35)adding  and  taking  the  square  root  of  the  sum,  i.e.. 


c,')2  P = (2(OS2(, 
C OS  0O 

tan2  0O  =C  2SINJC, 


and 

SIN  C , = SlN0o/(ctn2  P + SIN2  0O)1/2 
COSC,  = ctn/3/fctn1  P + SIN20O),/2 


I 


(D-37j 
(D- 38) 
fD-39) 


(D-40) 


To  proceed  with  the  analytical  solution  of  the  Lagrangian  equation,  the  time  dependency  of 
equation  (D-6)  /nust  be  eliminated.  Rewriting  the  equation  to  reflect  a conservative  system  (as 
previously  required)  results  in 


d7  - r02  - (r(OS2  0lX2  + **  = 0 
dt2  r2 

Introducing  the  dimensionless  expression  (ratio  of  the  earth’s  equatorial  radius  to  the  equatori  I 
component  of  the  missile's  spatial  position) 

Ws  R/rCOSfi  (D-4 1 ) 

the  first  and  second  time  derivatives  of  r,  (with  the  expression  for  the  angular  momentum)  becomes 

d 11  d / R \ 

dt  r2 COS2  0 dX  VWC'OS 0/  (D‘42^ 

and 

d2  r _ 11  d / dr  \ 

dt2  r2('()S20  dX  V dt  / 
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resulting  in 


d2  r 

H 
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/ H 

11  ( 

R 

dt2 

r2C  ()S30 

dX 

1 ( R/W)2 

dX  V 

W(  OS  0 

d2  r _ 

II3  W3  | 

W d3 

SIX  0 - 

SIX  0 

d3Wi 

dt2 

K3  \ 

dX3 

dX3  ) 

f equation  (D-9) 

may 

be  rewritten 

as : 

[amt)]’ 

Uy  substituting  the  value  for  dX/dt,  derived  in  equation(D-14)  ,the  above  equation  reduces  to 

r03=r(dOV  ( 1,2  ) 

VdX/  \ ( R/W)3COS0  / 


resulting  in 


r03 


n3  w3 


six  o 


(g) 


I'he  third  term  in  equation (D-9)  .reduces  to 


• , II2  W3 

(rCOSJ0)X2=  Rr  tosfl 

Substituting  the  above  results  into  equation  (D-9)  yields 

SEC  8 £*  -»  {at  " (g)’ 

O A 

- c os  o | = g cos3  o 

where  the  quantity  II 2 W2 / R3  was  factored  out  and  by  definition 

g = pR/H2 


(D-43) 


(U-44  ; 


(D-45) 


(D-46) 


(D-47) 


(D-48) 


Dividing  the  above  equation  (D-47),by  t le  SEC  0 and  evaluating  the  derivatives  of  the  SIX'  0, 
yields  the  following  linear  differential  equation  of  2nd  order,  i.e., 


where 


d3W 

d\5 


+ 


W = 


1 

SIX'3  0 


(11-49) 


SEC  0 s 1 1 + Cj  SIN3  (X  + C,)l  W2 


(D-30) 


I he  general  solution  to  the  above  non-homogeneous  linear  equation  is  made  up  of  the  sum  ot 
the  complementary  function  and  the  particular  integral.  1 he  solution  involves  two  independent 
arbitrary  constants  and  is  given  by 


W = Wc  + Wp 


(D- 50) 


D-7 


J 


where  W = We  is  the  general  solution  of  the  homogeneous  equation  (right-hand-side  of  equation 
(D-49 ) replaced  by  zero)  and  W = Wp  is  the  particular  integral.  Following  the  procedure  outlined 
in  any  differential  equations  textbook,  the  general  solution  is  given  by 

W=  (|  J |I+C2  SIN2  (X  + ( ,)|'/2+K,SIN(X-M  ,)  + 

K2  COS(X+C,  | 

where  K I . K2  are  constants  of  integrations  to  be  defined  and  C',  and  ( 2 are  the  integration 
constants,  previously  defined  by  equations  (D-33)  and  (D-27)  . 


(D-51) 


I o evaluate  the  undefined  integration  constants  K, , k2 , differentiate  equation  (0-51) , with 
respect  to  the  independent  variable  X which  is  the  longitudinal  coordinate  of  the  missile.  Multiplying 
W by  the  SIN  (X  + C, ) and  its  derivative  dW/dX  by  the  COS  (X  + C,  ),  adding  the  two  resulting 
equations  and  solving  for  k,  (by  rearranging  the  terms),  yields 

K,  = W SIN  (X  + C, ) + ^ COS  (X + C,)- 

cl  A 

0 + c2)  {•'  + c*  SINJ  <*  + <'.)l1/2  + 

SIN  (X  + C, ) + S|N  (XK'i  l COS^X+C, ) I 

II  +C2  SIN2  (X  + C,)|'/2  t CD-52) 

I he  k2  constant  of  integration  value  can  be  obtained  by  multiplying  equation  (D-51)  by  the  COS 
(X  + ( , ) and  dW/dX  by  sin  (X  + C, ).  Subtracting  the  two  resulting  equations  and  solving  for  k2. 
yields 

K2  = W COS  (X  + C, ) - dVV  SIN  tX  + C, ) - 

d CO-53) 

| SEC  0 COS  (X+C, ) - C?  SIN^X+C^)  COS  (X+(  , )l 
(l+Cj)  l 1 SEC  0 ( 

Imposing  the  previously  defined  set  of  initial  boundary  values,  i.e.. 


r r0,  6 -*  0o,  X -*  X0  -*  O 
the  above  constants  of  integration  reduce  to: 

dW’ 


and 


K,  = W()  SIN  C,  + 


dX 

dW 


COSC,  - g COS  0O  SIN  C, 


k2  - WD cost',  -£  SIN cI-(j^c  COS^COSC, 


(D-54) 
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Additional  mathematical  exercise,  enables  K;  to  be  written  as, 


K,  = Kj  tan  C,  + 


where 


and 


or 


, dW 

cg(-  r .?  

dX 

o ‘ (l+Cj) 

= R 

2 1/3 
|l  + C,  SIN  a2  1 ' 

dW 

dX 

dW 

dX 


‘o  ‘o 

_ r sec  e 


o 


COS  0 


_ W0r0 


° 1 tan  0o  COS  (3-  ctn  7| 
l v/cj*  SIN  Cj  COS  0 - ctn  y\ 


I o CDS  0 

The  analytical  solutions  to  the  set  of  2nd  order  ordinary  differential  equations  given  by 

-r02  -(rCOS2  0)X2  + 4 = C) 
dt2  r 

(r2  0)  + (r-2  SlN-2ii  X2  = O 

and  (r2  COS2  0)  X = constant 

can  now  be  summarized.  Since  the  system  is  conservative  (an  imposed  constraint),  equation 
(D-8)can  be  conveniently  integrated  to  obtain  the  total  extrinsic  angular  momentum  equation 

H = (r2  COS2  0)  X = constant 

H = r0  V0  COS  0O  SIN  y SIN  0 

that  is  functionally  dependent  upon  the  parameters  of  the  initial  state  vector. 


(D-55) 


f D- 56) 


(D-57) 


CD-9) 

CD-10) 

(D-8) 


(D- 11) 


With  the  above  relationship,  the  time  derivatives  in  equations  (D-9) , CO- 10)  can  be  eliminated, 
and  the  functional  dependency  of  0 on  X, yields  the  following  analytical  solution: 


or 

where 


0=0(X)  = arctan  | SIN  (X  + C,  )| 
0(X)  = ARC  SBC  ll+Cj  SIN  2 (X + (',)) 1/2 
C,  s arctar;  (SIN  0O  tan  0) 

/T?=  SEC  0O(  ctn  2 0 + SIN2  0O)‘/2 
( I + C, ) = C SC'2  0 SBC'2  0„ 


(D-29) 
CD- 30) 
(D-36) 
(D- 27) 
(D-58) 


D-9 


These  solutions  were  extensively  used  to  eliminate  the  time  derivatives  in  equation  (D-10) , 
thus  yielding  an  analytical  solution  that  relates  the  functional  dependency  ot  r on  X given  by; 


where 


V\  ( OS  0 = 


W 

II  +(  2 SIN2  (X  + ( , | 1/2 


and 


K, 

K2 


= K2  fan  C',  + 
= W0C’OSC,  - 


II  tCj  SIN2  (X  + ( ,)|1/2 


dW 

dX 


SECT, 

o 


1 

( I +C,) 


dW 

dX 


SIN  C , 
o 


g 

(l+Ca) 


+ K,  SIN  (X  + C , ) + 
KjCOStX  + C , ) 
SIN  C,  ( OS  0o 

COS  0o  COS  c, 


W0  = R/r0  C OS  0O 

dxL*(W()sS  I AT  SIN  C,  COSO-  ctnT| 

Note  that  the  above  solution  is  analogous  to  the  two-dimensional  ‘hit’  equation  derived  in 
Appendix  E. 


'D- 59) 

(0-51 } 

(D- 55) 

(D-54) 
(L)-  56) 
(D-57) 
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APPENDIX  E 

EVALUATION  OK  APPROXIMATE  VALUES  OF  THE  TWO  DIMENSIONAL 
BALLISTIC  MISSILE  PARAMETERS 


In  the  planar  representation  of  the  free-flight  trajectory,  the  equations,  derived  from  the 
Lagrangian  function,  are  defined  by  the  in-plane  polar  coordinates  (r,  0),i.e.,  A = O.  Thus,  for  a 
predetermined  surface  range,  constraints  are  imposed  on  the  burnout  parameters  ( v0.  ru  y)  to 
assure  this  objective.  In  the  following  derivations,  these  parameters  will  be  shown  to  be  functionally 
dependent  on  this  range.  The  in-plune  representation  of  equations  (D-9)  and  (D-10)  is  given  by: 


if*  e)-o 

dt 


(E-l) 


The  second  equation  expresses  the  conservation  of  angular  momentum  in  the  central  force  field 
and  may  be  expressed  as: 

H = i2  0 = r v SIN  7 3 rD  v0  SIN  y (E-2) 


where  the  burnout  velocity  vector  is  specified  by  its  speed  (vu)  and  flight  path  angle  with  respect 
to  the  local  vertical  (7).  Thus,  at  burnout  i.e.,  t = 0,  the  velocity  components  in  the  spherical  polar 
coordinate  system  (r,  0,  X),  at  r = rG,  0 = 0O  and  X = 0 are  given  by: 


dr 

dt 


= vG  COS  7 


d0 

dt 


= V°  SIN  7 
ro 


dX 

dt 


0 


(E-3) 


where  the  angles  are  illustrated  in  f igure  E-l . 


Since  the  solution  to  the  Lagrangian  equation  requires  finding  r and  0 as  functions  of  time, 
this  dependency  may  be  eliminated  by  rewriting  the  equation  of  motion  into  a differential  form 
i(0).  thus  eliminating  the  parameter  t.  f rom  the  expression  for  angular  momentum  (equation  (E-2) , 
the  differential  change  dt  is  related  to  a corresponding  change  d0: 

d = d0  d = H d 

dt  dt  dO  “ r1  dfl  (E-4) 


E-l 


and  the  second  derivative  with  respect  to  t and  0 is 


d2  _ H d /'ll  d \ 
dt2  r2  dO  \ r2  d 0/ 

The  Lagrangian  equation  becomes 

H d /Ji  dr  \ _ H2 

r2  dO  V r2  dO  / r3  r2 

Since 

1 =_d  n/r) 

r2  dO  dO 

and  by  introducing  a new  variable  u(0)  = l/r(0),  the  different!  i equation  for  the  orbit  becomes 


or 


H2  u2 


d2  u 
dO2 


+ H2  u3 


= - ptu 


2 


d2  u 
dO2 


+ u 


II2 


(E-5) 


(L-6) 


(E-7) 


(E-8) 


(E-9) 


with  a solution  given  by 

u=  K,  COS  ( K2  -01+m/H2 
du/dO  = K,  SIN  (K2  -0) 

where  K, . K2  are  constants  of  integration  solvable  with  the  initial  conditions 

uG  = u(O)  = 1 / r0 
11  o = " ro/ro2 


du 

= duo  - 

. ^o 

dO 

o 

o 

ii 

The  above  expressions  are  solved  for  the  constants  of  integrations. 

and  K,  SIN  K2  = - 

fo/H 

•••  \ = 

( 1 - M 

) COS  o - r° 

\ H2/ 

' ro  tl2 

) H 

which  reduces  to  the  expression 

u(0t  ru  = 

1 -COS  0 

+ SIN  (7  - 0 ) 

C SIN2  7 

SIN  7 

(E- 10} 


(E- 11) 


(F--12) 


(E- 13) 


F.-3 


1 he  parameter  is  a dimensionless  ratio  of  twice  the  initial  kinetic  energy  to  the  potential  energy 
at  burnout: 


l =v0J/<M/ r0i 


(E 


I he  burnout  speed  is  given  by: 

v ) = = M f 11  - C OS  0)  1 (E 

rQ  rQ  (_r0  u SIN2  7 + SIN  (0  - 7)  SIN  7J 

lo  insure  that  a ballistic  trajectory  will  actually  impact  at  a predetermined  range,  a particular 
combination  of  the  burnout  parameters  ( v"(>,  f0,  yi  is  necessary  to  achieve  a specific  range  angle 
0HR  . When  the  ballistic  missile  trajectory  intersects  the  earth's  surface  the  above  expression 
(equation  E - 1 3 ) becomes  known  at  the  “hit”  equation,  i.e.. 


Ir0!  _ I -COS0BR  + SIN(7~0BRI 
RfJ  I 1 SIN*  y SIN  7 


1 he  burnout  speed  required  to  impact  at  a range  0BR  , is  given  by. 


i?0i*=  « r 1 -COS0BR 

lfQl  |_(lr0l/lf,  I ) SIN2  7 + SIN  (0HK  -7)  SIN  7] 

where 


and 


= spherical  earth  radius  vector 


0|1K  = in-the-plane  range  angle, 

related  to  the  surface  range  S by 


(E 


(E 


0KR  = S/60.039  (DKG.) 

0HR  (DEG)  = lrj.:  I0BR  (RAD)/ 60. 039  fE 

Specifying  a set  of  burnout  velocities  (incremented  by  2000.  ft/s  in  magnitude  and  by  2.0  degrees 
in  direction)  at  a specific  burnout  altitude  ( I r(>  — r,  I ) the  free  flight  surface  range,  0 BR  . was 
calculated  and  plotted  as  a function  of  the  burnout  reentry  angle,  7.  Illustrated  by  Figure  E-2,  eacl 
iso-velocity  curve  ( 15,000.  ft/s  to  25,000  ft/s)  achieves  a maximum  surface  range  at  a specific- 
value  of  7 and  denotes  a minimum  energy  condition. 

lo  generate  these  minimum  energy  trajectories  requires  the  knowledge  of  the  magnitude  and 
direction  of  the  missile  burnout  velocity  vector  that  w ill  maximize  its  surface  range.  Given  the 
desired  surface  range  of  the  missile  0(tR  . the  functional  dependence  of  the  reentry  angle  7.  is 
obtained  by  setting  with  respect  to  7.  the  partial  derivative  of  equation  E- 16 , i.e., 

O0,jR/d7)  = O (i 


U-4 


f igure  E-2  Range  Angle  I Degrees)  versus  Burnout  Reentry  Angle  (Degrees)  for  a 
Family  of  Constant  Burnout  Velocity  Magnitude 


E-5 


resulting  in 


or 


tun  27  = SlNflHR/((  OS0hR  - lf0l/  I?,  I) 


2 7 = rr  - arctan 


sin  yltR 

1‘rj  I ) - COS  0 HR 


\ 


With  these  values  of  the  basic  parameters,  the  calculation  of  the  magnitude  of  the  burnout 
velocity  can  be  obtained  from  equation  E-20.  Thus,  fora  given  value  of  burnout  altitude,  these 
calculated  values  are  illustrated  in  Figure  E-3. 

To  clarify  and  bound  the  previously  introduced  dimensionless  parameter  C (defined  by 
equation (E-14)the  total  energy  of  the  conservative  system  is  given  by  the  sum  of  the  kinetic  and 
potential  energies,  i.e.. 


F = T + V = constant 


where:  T and  V are  the  total  extrinsic  kinetic  and  potential  energies,  respectively.  Since  the 
TOTAL.  energy  is  invariant  to  position,  the  above  equation  may  be  evaluated  at  the  apogee  and 
perigee  positions  in  the  interest  of  mathematical  expediency.  Therefore 

, _ 1 v a 1 p _ 1 V_T ' _ — 

2 ' i fa  I 2 V 

where  lv*al,  Inland  1 r*t  I . ifpl  are  the  magnitudes  of  the  velocities  and  positions  at  apogee  and 
perigee,  respectively.  From  the  general  expression  lor  the  total  extrinsic  angular  momentum  ot 
the  conservative  system,  7 x "v  at  apogee  and  perigee 


lral  KTal  = Ifpl  I v“p  i = constant 


resulting  in 


and  the  energy  equation  may  be  rewritten  as: 


IVpl 


{,,  S'l 

l I v „ / 


or 


CE-20) 

(E-21) 

CE-22) 

CE-23) 

(E-24) 

(K-2S) 


fb-2b) 


. . [ = 


2 


p (|  _ lvp|2|rpl 

Ifpl  \ 

! U - " 

irjji  \ i + rfpi/i?ai>  (irai  + irpi)  ) 

Thus,  the  total  energy  for  any  position  on  the  closed  elliptical  orbit  depends  only  upon  the 
magnitude  of  the  apogee-perigee  position  vectors,  since  the  total  energy  is  negative  and  since  the 
total  energy  at  burnout  may  be  written  as: 


(L-27) 


(E-28) 


< -M  =>C<2 
lrnl 


( t - 29 ) 


where  fv0l  and  Ir^l  represent  the  magnitudes  of  the  missile's  velocity  and  position  vectors 
respectively. 
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APPENDIX  t 

GRAPHICAL  SOLUTION  FOR  A QUICK  ASSESSMENT  OF  THE  DETECTION  PROBLEM 

To  aid  the  analyst  in  determining  the  visibility  of  an  object  (whether  it  be  a satellite  and/or 
missile)  in  the  space-time  domain,  a rapid  method  of  plotting  the  earth  track  (latitude  and  longitude) 
for  a specific  sensor(s)  location(s)  was  derived.  Using  a Mercator  Projection  Map  of  the  earth 
(illustrated  in  Figure  P-1,  the  spatial  object’s  ground  trace  was  plotted  on  the  HP  9862A  Plotter. 
Also,  the  volumetric  coverage  for  a specific  sensor  location  and  its  associated  visibility  parameters 
was  plotted.  The  contour  (illustrated  by  Figure  F-2)is  generated  from  the  projected  spatial  points 
on  the  circumference  of  the  base  for  a right-circular  cone.  This  cone,  with  apex  at  the  sensor’s 
location,  is  generated  by  specifying  a set  of  visibility  (constraint)  parameters,  i.e.,  elevation  angle 
and  slant  range.  (Ref . 5). 


To  plot  the  visibility  contours  on  the  niercator  projection  map  requires  a set  of  transformation 
equations  that  will  define  the  geographical  arc  length  between  the  set  of  projected  spatial  points 
and  the  specific  sensor’s  location.  Referring  to  Figure  F-3  from  the  law  of  cosines  equation  for 
a spherical  triangle,  the  cosine  of  the  arc  length  is  expressed  by 


COS  a = SIN  4>s  SIN  d>p  + COS  (AX)  COS  4>s 
COS  4>p. 


(F-l) 


where  AX  s \p  - \s  = longitudinal  separation  of  the  projected  point  and  sensor,  4>p,  <1>S  = latitude 
of  the  projected  point  and  sensor,  respectively. 


Since  the  arc  length,  a is  equivalent  to  the  subtended  earth  angle  (illustrated  in  Figure  F-4)  , 
its  functional  dependency  on  the  sensor’s  detection  range  and  elevation  angle  is  given  by: 


a = arc  tan 


Rs  COS  E ) 
R,,  + Rs  SIN  F I 


(F-2) 


where  the  sum  and  differences,  i.e.,  a ± 4>sdefines  the  latitude  of  the  most  northerly  and  most 
southerly  projected  spatial  points.  Thus,  the  spatial  object  is  said  to  be  visible  from  a ground 
tracking  sensor,  if  when  viewed  from  the  sensor,  it  lies  at  least  E degrees  above  the  true  horizon  of 
the  sensor  and  within  the  specified  detection  range.  A graphical  display  of  the  above  equation,  in 
I igure F-5  represents  a sensor’s  viewing  angle,  a (degrees)  as  a function  of  slant  range,  Ks  ( NM ) at 
detection  for  a specific  set  of  elevation  angles,  F (degrees). 


Since  the  Mercator  Projection  Map  can  only  display  a two-dimensional  volumetric  coverage, 
an  intersecting  ground  trace  (object’s  earth  track)  can  indicate  a possible  detection  provided  the 
minimum  altitude  requirement  is  met.  Referring  to  Figure  F-4, he  sensor’s  viewing  angle  may  be 
expressed  by; 

i R,  C OS  Ft 

a = arc COS^  R|  jrpj-J-  E (F-4) 
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Visibilit\  Contours  for  Specific  Sensors 


Figure  F-3  Three  Dimensional  Representation  of  the  Subtended  Earth  \ngle 
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Figure  F-4  Sensor  Volumetric  Coverage  Geometry 

and  is  plotted  (See  Figure  F-6)  as  a function  of  altitude.  H (NM)  for  elevation  angles  of'O..  :.0 
and  5.0  degrees.  I Ims,  the  detection  of  an  object  is  conditionally  dependent  upon  a specified  set 
(unprimed)  of  constraint  parameters.  I or  example,  it  an  object  is  confined  to  lie  within  a sphere- 
capped  cone,  its  parameter  v t (primed)  is  described  by; 

< «s  I 

F'  > E | (F-4) 

H'  5*  II 


F-5 


w 
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VIEWING  ANGLE,  a (DEG) 


Graphically , the  visibility  region’s  upper  limit  is  bounded  by  the  Rs  = constant,  curve  and  its 
intersection  with  the  constant  I curve  forms  the  lower  (constant  altitude)  boundary  value. 
Generally,  the  constraint  parameter  set  was  limited  to:  2.0  degrees  elevation  and  2000  NM 
detection  range  in  the  coverage  analysis  (See  f igure  F-7) . 

The  above  equations  were  programmed  for  the  III’  9820A  calculator  and  plotter  on  a 
mercator  projection  map  by  the  programmed  IIP  9862A  Plotter.  Referring  to  Figure  F-8  a general 
description  of  the  operational  sequence  is  illustrated  with  a detailed  description  published  in 
a RAIX' Technical  Memorandum.  (Ref.  5) 

The  computational  logic  follows: 

1.  Specify  the  desired  set  of  constraint  parameters  i.e.,  elevation  angle  (degrees)  and  slant 
range  (NM ) 

2.  Calculate  the  viewing  angle  a,  defined  by  equation  ( F-2)  . 

3.  Specify  the  sensor’s  location,  i.e.,  latitude  (0.  to  ± 90  degrees),  longitude  (20.  to  360. 
degrees,  I ).  For  example,  Xs  = 70°W  would  be  inputted  as  290°  I Sensor  location  less  than  20 
degrees  1 cannot  be  accommodated  on  this  particular  map  projection. 

4.  Calculate  the  latitude  of  the  most  southerly  projected  point,  i.e.,  <1>{,  = 4>s  - a 

5.  Calculate  the  longitude  of  the  projected  point,  defined  by  equation  ( F-l) . 

6.  The  sequence  is  repeated  (by  incrementing  the  latitude  of  the  projected  point,  lt>p  + A) 
until  the  visibility  contour  is  closed.  Thus  the  latitude  and  longitude  ot  a set  of  projected  points 

are  translated  into  the  X-Y  coordinate  system  by  the  Mercator  Projection  transformation  equations. 
These  equations  relate  the  linear  dependency  of  the  abscissa  on  longitude  and  the  non-linear 
dependency  of  the  ordinate  on  the  latitude  of  these  projected  points.  1 his  relationship  is  evident 
from  the  grid  scaling  illustrated  in  Figure  F-2. 

To  plot  the  earth  trace  of  a specific  trajectory  (an  output  of  the  Trajectory  Simulation 
Computer  Program)  the  same  set  of  transformation  equations,  relating  latitude  and  longitude  to 
t1,  ibscissa  and  ordinate  of  the  Mercator  Projection  Map.  were  programmed  and  extensively 
e,  oyed  in  the  analysis.  A detailed  description  of  the  program  and  equations  will  be  contained 
in  me  above  referenced  memorandum. 


Figure  F 7 Slant  Range  versus  Altitude  for  a Given  Value  of  Elevation  Range 
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